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Fingerprints of Dirac points in first-principles calculations of scanning tunneling 

spectra of graphene on a metal substrate 

J. Slawinska^'^ and I. Zasada^ 

^Department of Theoretical Physics and Computer Science, 

University of Lodz, Pomorska 149/153, 90-236 Lodz, Poland 
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Graphene physisorbed on a metal has its characteristic Dirac cones preserved in the band- 
structure, but the Fermi level of the system is shifted due to the interaction with the substrate. 
Based on density functional calculations with van der Waals corrections, we present a method to 
determine the position of the Dirac point with respect to the Fermi level from the measured scanning 
tunneling spectra (STS). It has been demonstrated that the dips in both simulated local density 
of states and in the observed dl /dV profiles are indeed the fingerprints of the Dirac points. The 
type and the level of doping can be then inferred directly from the STS data without any additional 
experimental technique. Test calculations of graphene on a Cu(lll) substrate have shown that the 
predicted position of the Dirac point is in close proximity to the experimental value reported in 
the recent studies. Moreover, simulations for graphene on a Pt(lll) surface allow us to explain the 
apparent contradictions in the state-of-the-art experimental works. 

PACS numbers: 73.22.Pr, 71. 15. Mb, 68.37.Ef 



I. INTRODUCTION 

The unique properties of graphene^, such as its ul- 
trathin geometry^ and high carrier mobilitj^ make it 
a promising candidate for applications in future nano- 
electronics, sensors and photonics^i^. Making complex 
devices requires the production of large enough high- 
quahty graphene sheets to achieve scalability and ben- 
efit from exceptional electronic properties observed in 
the flat domains of the monolayer sheets^. In practice, 
the large-area graphene has been grown epitaxially on 
transition metals ^"^ , especially the chemical vapor depo- 
sition (CVD) method has been developed to synthesize 
graphene on copper— and gold^^ substrates, although 
their catalyst's role is still not well understoodi^. Due to 
the high quality and transferability of the samples pre- 
pared on metallic surfacesi^, much theoreticali^"— and 
experimental eSoi^^'^— has been devoted to shed light 
on the mechanisms of graphene- metal interactions as well 
as on the modifications of its electronic properties. In 
particular, the doping effect^^ in physisorbed graphene 
has been widely studied both theoretically and experi- 
mentally by using a number of techniques, as for example 
density functional theory (DFT) , scanning tunneling mi- 
croscopy (STM) and angle-resolved photoemission spec- 
troscopy (ARPES). 

Scanning probe methods^^ have developed into a pow- 
erful tool for the determination of the structure of sur- 
faces and interfaces, thus it is especially suitable for stud- 
ies of graphene deposited on conductive substrates. The 
STS mode allows us to record the first derivative of the 
tunneling current with respect to the bias voltage which 
is the measure of the surface density of states at every 
arbitrarily fixed point. The electronic structure is probed 
locally^^, thus the presence of heterogeneity at different 
places on the surface or periodically repeated domains 
can be identified. Moreover, the information as to how 



the metal surface states are modified by the adsorption 
of graphene sheet can be obtained after careful analysis. 

According to early local density approximation (LDA)- 
based DFT studieaiS, of a graphene/Cu(lll) system, the 
physisorption leads to the change in the location of Dirac 
point, since the electrons are donated by a substrate to 
graphene (n-type doping of -0.17 eV). In a recent experi- 
mental study of graphene epitaxially grown on a Cu(lll) 
surface the results of STS measurements are reported^!. 
The dip in the presented STS profile (Fig.SC in RefHH) 
at approximately -0.35 eV is associated with the Dirac 
point shifted below the Fermi level, which qualitatively 
confirms previous theoretical results. However, the inter- 
pretation of differential tunneling spectra is not straight- 
forward due to the influence of the tip states2I. Thus, 
as suggested in Refl2ll only a combination of STS with 
a secondary experimental technique could give a direct 
evidence that a dip reflects the local density of states 
at the Dirac point of graphene. ARPES measurement 
could conflrm the energy position of the Dirac point, 
however sufficiently large monocrystalline samples are re- 
quired to find specific paths inside the Brillouin zone of 
the system. On the other hand, the new methods of 
Dirac point mapping have been proposed, such as fol- 
lowing the gate-dependent position of adjacent dips in 
the tunable system??^ or analysis of Fourier transforms of 
STM maps revealing a linear dispersion relation for states 
above and below the Fermi energy- -. They are, however, 
either not applicable for conductive substrate or at least 
hardly achievable without special equipement such as low 
temperature (LT) STM with lock-in detection. 

In contrast, the theoretical predictions and the exper- 
imental findings for the graphene/Pt(lll) system seem 
to be in an excellent agreement. The adsorption distance 
of about 3.30 A is both predicted by DFT calculationgi^ 
and measured by low-energy electron diffraction (LEED) 
technique.^ The doping level is estimated to approxi- 



mately 0.33 eV by DFT calculationgi^ while ARPES mea- 
surements gave the value of 0.3 eV.^ However, ARPES 
can determine only occupied states: in the case of the 
Fermi level shift downward with respect to the Dirac 
point, the value of doping must be extrapolated. In 
this view, the main advantage of STS is its ability to 
record information about unoccupied part of the spec- 
trum, which allows us to measure doping directly. Recent 
experimental studies revealed two types of spectra for 
the graphene/Pt(lll) system: the one with clearly visi- 
ble onset of a Pt (111) surface state^i and the second one 
with global minimum corresponding to the Dirac point^ 
(in agreement with calculated PDOS). We propose the 
possible explanation of this difference due to analysis of 
metallic surface states in terms of the LDOS above the 
surface. 

All these indicate that linking theoretical predictions 
and experimental findings needs an approach indepen- 
dent of the limitations of measurement and preparation 
methods. In this paper, we report on the first-principles 
simulations of differential spectra that can be directly 
compared with measured dl /dV profiles of graphene on 
a metal surface. As the DFT method provides the infor- 
mation about the band structure, the dip in calculated 
LDOS above the surface can be unambiguously linked 
to the energy position of Dirac points in the dispersion 
relation. The examples of graphene on a Cu(lll) and 
Pt(lll) substrates have been carefully elaborated to test 
the presented method. 

The paper is organized as follows: details of DFT cal- 
culations are described in Sec. II, whereas the idea of a 
proposed method as well as two examples of its applica- 
tion are presented in Sec. III. Final remarks are given in 
Sec. IV. 



II. DFT CALCULATIONS 

All DFT calculations have been performed using VASP 
package^Si^ equipped with the projector augmented 
wave (PAW) method^^'^'^ for electron-ion interactions. 
The exchange-correlation energy is calculated using the 
generalized gradient approximation within the Perdew, 
Burke and Ernzerhof (PBE)^^ parametrization scheme. 
Long-range dispersion corrections have been taken into 
account within a DFT-D2 approach of Grimme^, as im- 
plemented in the latest version of VASP"^-. The pair in- 
teractions up to a radius of 12 A have been included in 
the calculations and the global scaling factor sg bas been 
set to 0.75 due to the choice of the PBE functional. The 
dispersion coefficients Cg and van der Waals radii Rq for 
C and Cu atoms are defined in the code, according to 
the suggestion of Grimme^, while the Cg and Rq coef- 
ficients of Pt have been set to 20 Jnm^/mol and 1.9A, 
respectively, as in Ref. l20li^ The electronic wave func- 
tions have been expanded in a plane-wave basis set of 
400 eV, while the electronic self-consistency criterion has 
been set to 10~^ eV. 
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FIG. 1: (Color online) Top view (a) and side view (b) of ad- 
sorption geometry of graphene on a Cu(lll) substrate. Car- 
bon atoms are denoted as red (darker) balls, copper atoms 
as gray balls. Top view (c) and side view (d) of adsorption 
geometry of graphene on a Pt(lll) surface. Platinum atoms 
are denoted as yellow (lighter) balls. Parallelograms define 
the unit cells. 



The (111) surfaces of metals are simulated by a peri- 
odic slab geometry. Supercells containing six (five) lay- 
ers of Cu (Pt) atoms and a graphene sheet adsorbed on 
top are separated by the vacuum thickness of at least 
20 (15) A. The typical in-plane adsorption geometry is 
shown in Figlli;a),(b) and Figd^c) , (d) for Cu(lll) and 
Pt(lll) substrate, respectively. The graphene-metal in- 
terface is modeled by direct matching of a graphene's 
and a substrate's unit cells, although in the real sys- 
tems the lattice-mismatch leads to the Moire pattern on 
the surface. This artificially generated strain tends to 
affect the surface relaxation and the electronic proper- 
ties, which were discussed in Ref. [l^. Thus, for the 
graphene/Cu(lll) system we have performed calcula- 
tions for both configurations: with the in-plane lattice 
parameter ahex equal to the PBE+D2-optimized value 
for graphene 2.466 A as well as with the lattice constant 
adapted to that of a metal (two cases were considered: 
PBE+D2 optimized ahcx =2.524 A and an experimental 
value of 2.56 A). There exist also three inequivalent stack- 
ing orders of graphene and Cu layers. The carbon atoms 
can lie above metal atoms in layers 1 and 3 (known as 
the top- fee configuration, see FigH]), 1 and 2 (a top-hep) 
or 2 and 3 (a hep- fee). The previous LDA studyiS iden- 
tified the top-fee configuration as the most stable one, 
which was also confirmed by the present vdW calcula- 
tions (total energies of top- hep and hep- fee configurations 
are higher than top-fcc one by 2.3 meV and 13.1 meV, 
respectively). In view of the above test calculations, in 



the case of graphene/Pt(lll) system we have used only 
one value of in-plane lattice paramater ahcx a-s well as 
2x2 graphene's unit cell in configuration defined in Fig. 
[IJc). During all structure relaxations the metal in top 
two layers as well as all carbon atoms have been allowed 
to move. Total energies were converged to within fO^® 
with respect to the ionic steps. 

The problem of accurate Brillouin zone sampling is of 
particular importance for simulations of local density of 
states above the surface. The simulations of STM images 
are not computationally demanding, but the resolution 
of the calculated STS profiles is determined by the accu- 
racy of the band-structure maps. For example, the iden- 
tification of the onset of a noble metal surface requires 
at least one k-point every few meV for a specific band 
which means a few thousands of k-points in the Brillouin 
zone^^. During the relaxation and the self-consistent run 
we applied the tetrahedron scheme^^ and the F-centered 
36x36xf and 24x24x1 k-point meshes for small and 
large supercell, respectively, whereas in the final local 
density of states calculations, a dense 336 x 336 x 1 
(192 X 192 x 1 for larger supercell) F-centered k-points 
mesh and the Gaussian smearing of width a =0.05 eV 
have been used. 



III. IDENTIFICATION OF DIRAC POINTS IN 
STS DATA 

The Tersoff-Hamann model of tunneling, due to its 
simplicity and qualitative reliability, is implemented in 
VASP as well as nearly every other DFT code to make 
simulations of pseudo-STM topography images. In this 
approach the tunneling current is proportional to the lo- 
cal density of states at the position of the STM tip^: 
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Provided that the k-points mesh is sufficiently fine, the 
reliable pseudo-STS profiles can be obtained by an eval- 
uation of the constant-height charge images for several 
values of bias voltage followed by their numerical differ- 
entiation. The local density of states is automatically 
calculated in every single point above the surface and 
the choice of a spectrum in the particular position could 
correspond to the STS data taken under open loop con- 
ditions at the fixed point. 



A. Cu(lll) substrate: n-type doping of graphene 

The local densities of states of a clean perfect Cu(lll) 
surface are shown in Fig. [31 We have evaluated the spec- 
tra for all three considered choices of the in-plane lattice 
paramater ahox- Their overall shape is similar to those 
presented in the previous theoretical studie a^^i^° as well 
as to the data obtained in a number of experiments^^—. 



The step-shaped onsets of the surface state can be easily 
observed at bias voltages Ui — —700 mV, U2 — —440 
mV and [/a — —280 mV for the values of lattice con- 
stants a^p^=2.466, a^(,^=2.524 and a^j,^=2.56 A, respec- 
tively. It should be noted that the oscillations present 
in the profiles originate from a discrete k-points mesh. 
According to the discussion presented in Ref.^^, the sig- 
nature of a surface electronic structure should be very 
sensitive to the choice of the lattice constant. Indeed, 
the experimental position of the surface state is repro- 
duced only for a value of the lattice constant a^^^ opti- 
mized within the PBE+D2 approach. This appears to 
be in contrast with previous theoretical results based on 
a Perdew-Wang (PW91) parametrization, where setting 
the ahcx to an experimental value was needed to obtain 
the quantitatively correct data. The discrepancies seem 
to be caused by using of a different parametrization ap- 
proach and will be widely discussed in a separate study. 
The problem of the surface onset shift is especially im- 
portant for modeling the graphene-metal interface where 
adjusting the substrate's lattice constant to graphene is 
usually requiredi^. It means that a correct theoretical 
determination of subtle features such as experimentally 
observed suppression of a Cu(lll) surface state in the 
presence of graphene^i could be hardly achievable within 
a 1 X 1 supercell approach. 

The bonding of graphene to a copper surface is rather 
weak and the charactersitic Dirac cones are nearly pre- 
served in the band-structure. Both LDA'*'*^ and the 
present DFT-D2 calculations predict a tiny band-gap of 
the order of 10 meV in the case of the matched unit cells, 
but in a real incommensurate interface it should disap- 
pear. It is worthwile to notice the interesting fact that the 
presence or lack of this gap is closely related to the sym- 
metries (hexagonal or trigonal) being preserved in the 
system. Dirac fermions can exist even in chemisorbed 
graphene provided that the defects are distributed sym- 
metrically. This problem is widely and elegantly dis- 
cussed in Ref. |45| in terms of symmetries, group theory 
and tight binding Hamiltonians. It has been also shown 
by means of DFT calculations^ that despite functional- 
ization through dopants or vacancy defects, the presence 
of Dirac fermions may be recovered for superstructures 
and nanomeshes with specific symmetry. In the case of 
Cu(lll) substrate, however, real weak interaction with 
graphene has been confirmed experimentally. 

The physisorption of graphene on a Cu(lll) sub- 
strate causes a Fermi level shift upward--. In 
FigE] (a) we present the energy dispersion relation of 
graphene/Cu(lll) system for a value of ahcx =2.466 A. 
The similar band-structures have been obtained for other 
values of the lattice constant, i.e. 2.524 A and 2.56 
A. The optimized adsorption distances are about 3.06 
A 2.98 A and 2.88 A, respectively. It is worthwhile 
to note that DFT-D2 tends to overestimate binding en- 
ergies, e.g. for the considered system binding energy 
per carbon atom equals 117 meV, which is very close to 
the values for physisorbed graphene reported in a recent 
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FIG. 2: (Color online) The simulated LDOS spectra of a clean 
Cu(lll) surface for different choices of the in-plane lattice 
parameter ahcx- 



DFT-D2 study—, but significantly differs from 38 meV 
calculated within vdw-DF approach by Vanin et al^. 
Next, we define the doping level as AEf= Ed — Ep, 
i.e. a difference between Dirac point energy and a Fermi 
level. The values of doping levels are equal to -0.427 eV, 
-0.584 eV and -0.686 eV for the considered lattice param- 
eters, respectively. The problem of changing in doping 
with the variation of the in-plane lattice parameter has 
been discussed in details in Ref. 19 : it seems that the 
most reliable value of doping is provided by the choice 
of in-plane lattice parameter matched to graphene's one. 
The differences in doping can be also described in terms 
of the charge transfer, in particular evaluated using the 
Bader scheme^"—. The charge transfer SQ equals to - 
7.4, -24.1 and -31.3 xl0~^e, respectively, where negative 
sign indicates electron transfer to graphene. The increase 
in SQ with larger values of ahcx seems to be related to 
decrease in adsorption distance. Moreover, the adsorp- 
tion distances are usually underestimated^" by DFT-D2, 
thus the calculated doping level can be slightly below the 
measured values. 

The STM experiment of graphene on a Cu(lll) sub- 
strate revealed a Moire pattern present on the surface^l. 
It means that the length of graphene bonds is preserved, 
which seems to be common for all weakly binding sys- 
tems. Moreover, it has been previously demonstrated for 
a similar system, that the most reliable values of doping 
are obtained when a unit cell is adjusted to graphene's 
lattice constant. For example, in case of graphene - gold 
interface, stretching of bonds leads to the value of doping 
which is in disagreement with the experimental data^^. 
The choice of ahox= 2.466 A should then lead to the most 
accurate description of the graphene's electronic proper- 
ties. 
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FIG. 3: (Color online) (a) Left-hand panel: The band- 
structure of freestanding graphene. The inset shows a Bril- 
louin zone with the F, K, and M high-symmetry points. 
Right-hand panel: The electronic structure of graphene on 
Cu(lll). (b) The calculated LDOS profile above the surface 
of a graphene/Cu(lll) system for a sample-tip distance larger 
than 2 A (solid line) and equal to about 1 A (dashed line). 
The parameter Qhcx has been set to 2.466 A. The zero energy 
is at a Fermi level, (c) Simulated constant-height STM im- 
age corresponding to calculated LDOS presented above. Bias 
voltage is set to U =-100 mV. Light spots are in the positions 
of carbon atoms. 



In Fig. [3] (b) we present a LDOS spectrum 
for a graphene/Cu(lll) interface simulated with the 
ahGx=2.466 A (sohd hne). Since the position of the 
Dirac point is known from the band-structure (see Fig. 
[3] (a)), we can directly map it onto an STS profile. It 
means that a dip at -0.427 eV in the spectrum (Fig. [3] 
(b), solid line) can be associated with the Dirac point in 
the energy dispersion, which provides an unambiguous 
link between experimental STS profiles and the calcu- 
lated band-structures. In this light, the local minimum 
observed at -0.35 eV in the spectrum (Fig.SC in Rcf.^^) 
should be indeed interpreted as a signature of a Dirac 
point. The measured value of doping is quite close to the 
one predicted in this study. The difference of about 80 
meV might be explained by the effect of phonon-mediated 
inelastic tunneling as recently discussed in Ref.^^. This 
phenomenon could cause a systematic shift of a Dirac 
point dip toward the Fermi level by about few tens of 
meV. Another reason for a doping level shift may be con- 
nected with the adsorption distance underestimated in 
the DFT-D2 approach. 

Despite the oscillations related to the finite resolution 
of the spectrum presented in Fig. [31(b), a Dirac point dip 
can be still easily recognized. However, if the profile is 
evaluated using the minimal required k-points mesh the 
Dirac point can be hardly identified. One may addition- 
ally confirm its location by decreasing the tip-sample dis- 
tance, i.e. approaching to the surface (Fig. [21(b), dashed 
line) . At the distances of the order of 1 A above the sur- 
face, the graphene states dominate in the spectrum and 
the minimum can be unambiguously recognized. 

It should be also noted that the onset of a Cu(lll) 
surface state after the deposition of graphene (Fig(31(b), 
solid line) is shifted toward the Fermi level with respect 
to the position predicted in the calculations for a clean 
surface done using the same lattice-constant parameter 
ahcx=2.466 A (Fig. [21 , solid line). It appears to be 
connected with the changes in the work function upon 
adsorption of the carbon atoms. 



B. Pt(lll) substrate: p-type doping of graphene 

In the case of the Ft (111) substrate, we have restricted 
our analysis to only one fixed lattice constant (see RefUOl 
for wider discussion). In Fig. 4(a) we present the band- 
structure of graphene interacting with Ft (111): we can 
observe that, in accordance with previous result a^^'^" , 
Dirac cones are preserved and shifted upward with re- 
spect to the Fermi level, which results in p-type dop- 
ing of graphene (0.41 eV). In spite of the fact that we 
used only 2x2 unit cell of graphene the value of doping 
agrees with that presented in Refl2(]l where the Moire 
pattern was taken into account. It should be noted that 
the band-structure measured by ARPES(|^) also confirms 
the presence of Dirac fermions in this lattice-mismatched 
system. 

In Fig. 4(b) we present a LDOS profiles of the 
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FIG. 4: (Color online) (a) Left-hand panel: The band- 
structure of freestanding graphene in 2x2 unit cell. The 
inset shows a Brillouin zone with the F, K, and M high- 
symmetry points. Right-hand panel: The electronic struc- 
ture of graphene on Pt(lll). (b) The calculated LDOS profile 
above the surface of a graphene/Pt(lll) system for a sample- 
tip distance larger than 5 A (solid line) and equal to about 
3 A (dashed line) . (c) The simulated LDOS profile of a clean 
Pt(lll) substrate. It agrees well with previous theoretical 
and experimental spectrai^. The zero energy is at a Fermi 
level. 



graphene/Pt(lll) system and, for comparison, spectrum 
of the pristine Pt(lll) surface (Fig. 4 (c)). Similarly 
as in the case of graphene/Cu(lll) system, we can as- 
sociate the position of the Dirac point in the band- 
structure Ejj with a dip at energy 0.41 eV in calculated 
LDOS spectra above the surface. Linking two types of 
data allows to recognize the fingerprints of Dirac points 
in simulated and measured dl/dV spectra. However, 
one can notice at least two differences comparing with 
the graphene/Cu(lll) system: (i) the surface states of 
Pt(lll) appear to be better pronounced: even in the 
presence of global minimum responsible for Dirac point 
(dashed line in Fig. 4(b)) the surface state of platinum 
can be easily recognized; (ii) it is not unambiguous as to 
which of two dips in the LDOS curve (solid line in Fig. 
4(b)) reflects the Dirac point: one of them may be in- 
duced by Pt surface state; in this context the evaluation 
of LDOS close to the surface is needed (dashed line). 

In this light, it is not surprising that contradictory ex- 
perimental results have been recently reported. In STS 
data presented in Refl24l (Fig. IB therein) there is no 
signature that could be associated with the presence of 
Dirac point in any particular position. It is supposed^! 
that this is due to Pt(lll) surface state features which 
are still seen in the spectra. We suggest that in such cases 
only direct and exact comparison with simulated LDOS 
might allow one to identify the Dirac point. On the other 
hand, dl/dV spectrum shown in Refl20l (Fig. 1(e) therein) 
corresponds rather to the dashed line in Fig 4 (b) in this 
present paper as well as to PDOS profile calculated in 
Refl2C)l. In this case the surface states are hardly visible 
which allows one to associate the global minimum to a 
Dirac point. The absence or presence of surface state 
may be related to preparation technique or details of the 
STS experiment. 



IV. FINAL REMARKS 

In summary, we have performed DFT-D2 calculations 
for graphene interacting with Cu(lll) and Pt(lll) sub- 



strates and evaluated the corresponding LDOS spec- 
tra above the surface. The fingerprints of Dirac points 
present in the band structure are easily recognized in the 
simulated STS profiles, which enables a straightforward 
comparison of theoretical and experimental data. The 
approach can be successfully used even for the systems 
where the Dirac point and the surface state onset seem 
to be overlapped. The doping of graphene on Cu(lll) 
predicted in this study is in a better agreement with the 
measured value than the previous LDA results. It seems 
that using the dispersion corrections in DFT calcula- 
tions is essential for accurate determination of electronic 
characteristics and demonstrates that vdW forces are re- 
sponsible for properties of the graphene-metal interfaces. 
Moreover, the simulations for the graphene/Pt(lll) sys- 
tem explain the difficulties with identification of Dirac 
points in the experimental STS spectra. 
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